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PULSAR PARALLAXES AT 5 GHZ WITH THE VERY LONG BASELINE ARRAY 
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ABSTRACT 

We present the first pulsar parallaxes measured with phase-referenced pulsar VLBI observations at 
5 GHz. Due to the steep spectra of pulsars, previous astrometric measurements have been at lower 
frequencies. However, the strongest pulsars can be observed at 5 GHz, offering the benefit of lower 
combined ionospheric and tropospheric phase errors, which usually limit VLBI astrometric accuracy. 
The pulsars B0329-I-54, B0355-I-54 and B1929-I-10 were observed for 7 epochs spread evenly over 2 years. 
For B0329-I-54, large systematic errors lead to only an upper limit on the parallax (tt < 1.5 mas). 
A new proper motion and parallax were measured for B0355-I-54 (tt = 0.91 ± 0.16 mas), implying a 
distance of 1.04lQ^g kpc and a transverse velocity of 6llg^ km s~^. The parallax and proper motion 
for B1929-I-10 were significantly improved (tt — 2.77 ± 0.07 mas), yielding a distance of 36l!tg° pc and 
a transverse velocity of 177lg km s~^. We demonstrate that the astrometric errors are correlated with 
the angular separation between the phase reference calibrator and the target source, with significantly 
lower errors at 5 GHz compared to 1.6 GHz. Finally, based on our new distance determinations for 
61929+10 and B0355+54, we derive or constrain the luminosities of each pulsar at high energies. We 
show that, for thermal emission models, the emitting area for X-rays from PSR B1929-I-10 is roughly 
consistent with the canonical size for a heated polar cap, and that the conversion of spin-down power 
to gamma-ray luminosity in B0355-I-54 must be low. The new proper motion for B1929-I-10 also implies 
that its progenitor is unlikely to have been the binary companion of the runaway 0-star C, Ophiuchi. 

Subject headings: astrometry — pulsars individual (PSR B0355-f 54, PSR B1929+10) — 
stars:kinematics — stars:neutron — X-rays:stars 



1. INTRODUCTION 

The measurement of distances is a fundamental prob- 
lem in astronomy. Trigonometric parallax offers one of 
the only model-independent mechanisms for such measure- 
ments, thus defining the lowest rungs of the distance lad- 
der. The parallaxes and proper motions of pulsars, ob- 
tained through imaging or pulse-timing observations, pro- 
vide measurements of their distances and velocities, and 
thus allow us to probe the physics of neutron star forma- 
tion in supernovae, as well as the intervening interstellar 
medium which scatters and disperses the radio pulses from 
pulsars. 

Pulsar distances provide crucial calibration for electron 
density distribution models of the Galaxy (e.g. Taylor & 
Cordes 1993; Cordes & Lazio 2002), which, in turn, pro- 
vide estimates of the distances for the majority of the pul- 
sar population and underlie estimates of their velocities. 
The velocity distribution of pulsars (Arzoumanian, Cher- 
noff & Cordes 2002) helps constrain core collapse processes 
in supernovae (Burrows & Hayes 1996), while velocities of 
individual objects can be used to trace their progenitor 
supernovae and birth sites (Hoogerwerf, de Bruijne & de 
Zeeuw 2001). Recent results from X-ray observations, cou- 
pled with distance estimates, have established constraints 
on neutron star cooling curves (Tsuruta et al. 2002) which 
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limit plausible nuclear Equations of State. 

While the diverse applications of pulsar distances and 
velocities make a compelling case for the necessity of pul- 
sar astrometry, reliable measurements have been few and 
far between (see, e.g. Toscano et al. 1999). The paucity of 
parallax measurements is due to the difficulty of measuring 
the effect. For a pulsar at 1 kpc, the angular displacement 
is only 1 mas and the semi-annual timing perturbation is 
1.6 /iS. Few pulsars allow sufficient precision to measure 
the timing effect, while astrometry with ad hoc very long 
baseline interferometry (VLBI) arrays has been difficult 
and inconsistent. However, the advent of the NRAO Very 
Long Baseline Array, a full-time, dedicated VLBI instru- 
ment, coupled with pulsar gating and refined calibration 
techniques, has significantly altered the situation and led 
to several parallax measurements in recent years (Brisken 
et al. 2000; Chatterjee et al. 2001; Brisken et al. 2002). 
Due to the steep spectra of pulsars [Sy oc i^", (a) ~ —1.6; 
Lorimer et al. 1995) these measurements have all been at 
1.4 GHz, a frequency lower than that used for most other 
VLBI astrometry. Here we present the first pulsar par- 
allaxes measured at 5 GHz. In §2, we outline the ob- 
servations, calibration process and data analysis. Astro- 
metric results for the pulsars B0329-I-54, B0355-I-54 and 
B1929-I-10 are presented in §3, and compared to those of 
Brisken et al. (2002) for B0329+54 and B1929-H10. We es- 
timate the attainable astrometric accuracy at 1.4 GHz and 
5 GHz in §4, and discuss the implications of the measured 
distances and velocities for B1929+10 and B0355-t-54 in 
§5. 

2. OBSERVATIONS AND DATA ANALYSIS 

Astrometric VLBI observations require the visibility 
phase of the target source to be calibrated using obser- 
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vatons of a nearby calibrator source. Typically, the tar- 
get and the "nodding" calibrator are observed in short 
alternate scans. To achieve phase connection across mul- 
tiple scans, phase referencing requires rapid cycle times 
and a small angular separation between the target and the 
nodding calibrator (hereafter referred to as the "calibrator 
throw"). The limiting source of astrometric error is typi- 
cally the uncalibrated differential atmosphere between the 
target and the calibrator, which includes an ionospheric 
dispersive phase error (oc v~^) and a tropospheric non- 
dispersive phase error (cx u). These two error terms arc 
comparable in size at the 5 GHz band and provide an 
overall approximate minimum error, making the band well 
suited for astrometry. For example, a 0.36 mas parallax 
has been measured for Sco X-1 by Bradshaw, Fomalont, 
& Geldzahler (1999) at 5 GHz. We chose three pulsars, 
B0329+54, B0355+54 and B1929+10, which have flux 
densities between 1 and 10 mJy at 5 GHz, and observed 
them for 7 epochs separated by 4 months each between 
2000 March and 2002 March. Solar maximum occurred 
within this period, leading to heightened ionospheric ac- 
tivity and ~ 3-10 x increases in the daily average total 
electron content compared to that at solar minimum as 
measured by GPS (e.g. Warnant & Pottiaux 2000). Such 
ionospheric activity is expected to cause epoch-dependent 
systematic errors in our astrometry. 

At each epoch, the pulsars were phase-referenced to ex- 
tragalactic sources included in the International Celestial 
Reference Frame. B0329-I-54 and B0355+54 were phase- 
referenced to ICRF J035929. 7+505750 (NRAO 150), from 
which they are separated by 5.3° and 3.3°, respectively. 
For B1929-I-10, ICRF J192840.8-I-084848, a calibrator 2.3° 
away, was used as a reference source, and each obser- 
vation included a few scans on ICRF J193435.0-M04340 
(hereafter J1934-I-1043), only 0.6° from the pulsar. The 
absolute positions of the calibrator sources are known to 
between 2-3 mas (Ma et al. 1998). The observations alter- 
nated between 90 seconds on the calibrator and 150 sec- 
onds on the target, and the pulsar gate was applied to the 
data at the VLB A correlator. Gating boosted the S/N ra- 
tio by a factor of - 3-4, where / = To„/(To„ + Toff ) 
is the gate duty cycle. AIPS, the Astronomical Image Pro- 
cessing System, was used for the data reduction, which 
followed standard procedures (Beasley & Conway 1995). 
The basic steps include flagging of bad records, amplitude 
calibration based on the system temperatures at each an- 
tenna, fringe fitting, and several passes of self calibration 
on the nodding calibrator. The multiple self calibration 
passes serve to correct the visibility phases and gains while 
fixing the position of the nodding calibrator. 

For B0329+54 and B0355+54, the cahbration solutions 
derived for the nodding calibrator were transferred to the 
pulsar data, images were created, and the AIPS task JMFIT 
was used to obtain astrometric positions at each epoch. As 
described by Chatterjee et al. (2001), position uncertain- 
ties were derived as a combination of the synthesized beam 
resolution and systematic errors in quadrature, where the 
systematic errors were estimated at each epoch using the 
deviation of the deconvolved image from that expected 
for a point source. While the typical uncertainties are 
< 1 mas, the systematic errors were worse for B0329+54 
at every epoch, and exceeded 5 mas at one epoch (2000 



November, due to ionospheric activity) for both pulsars, 
rendering that epoch useless for astrometry. 

For B1929-I-10, after transferring the calibration so- 
lutions from the nodding calibrator, images were cre- 
ated for both the pulsar and the second reference source 
J1934-I-1043, and JMFIT was used as above. Ideally, a con- 
stant position should be obtained for J1934-I-1043 at each 
epoch, but the observed position varied randomly, with 
a median angular scatter (i.e., the median of the angu- 
lar position offsets) of 0.7 mas over 7 epochs. Phase self- 
calibration of J1934+1043 fixes its position at the expected 
location, and transferring these solutions to B1929-I-10 
not only shifts the observed pulsar position, but also im- 
proves the image quality and reduces the estimated sys- 
tematic errors significantly. Thus the astrometric posi- 
tions for B1929-I-10 were referenced to J1934-hl043, 0.6° 
away. Here, too, position uncertainties at one epoch (2001 
March) exceeded 5 mas, and we discarded it in the final 
astrometric fits. 

3. PARALLAXES AND PROPER MOTIONS 

The derived positions and position uncertainties for each 
pulsar at 5 GHz were used to fit for an annual trigono- 
metric parallax, along with zero-point position offsets and 
proper motions in each coordinate (5 parameters), us- 
ing a weighted least-squares analysis. For B0329+54 and 
B1929-I-10, positions measured by Brisken et al. (2002) at 
1.4 GHz were also used to perform joint fits: two inde- 
pendent zero-point position offsets were employed for the 
two data sets to account for the different nodding calibra- 
tors and observing frequencies employed. We describe the 
results for each pulsar in turn. 

For B0329-I-54, unmodeled systematic errors dominate 
the fit and no satisfactory solution was obtained. The for- 
mal value of the parallax is tt = —0.45 ± 0.36 mas, with 
a reduced = 2.1 for 7 degrees of freedom. Holding the 
parallax fixed at does not improve the fit for proper mo- 
tion, yielding Ha = 16.8 ± 0.4 mas yr~^, ns = —11.1 ± 
0.9 mas yr"-'^, and a reduced = 2.3 (8 degrees of free- 
dom). Given the position uncertainties, we estimate an up- 
per limit on the parallax tt < 1.5 mas, which is consistent 
with the Brisken et al. (2002) value of 0.94±0.11 mas. In a 
joint fit to both data sets, the larger position uncertainties 
for the 5 GHz data leave the earlier parallax unimproved. 
The joint proper motion is also within the errors of the 
earlier Brisken et al. value of /i^ ~ 17.00 ± 0.27 mas yr~^ 
and us = -9.48 ± 0.37 mas yr"!. The 5 GHz data do 
not contribute useful astrometry in this case, because the 
calibrator was too far away (5.3°) to adequately correct 
phase perturbations in the data, as we discuss below. 

The position uncertainties are better behaved for 
B0355-I-54 and a good weighted least squares fit was ob- 
tained for the parallax and proper motion, with a reduced 
= 0.82 for 7 degrees of freedom. The derived best-fit 
values are listed in Table 1, and imply most probable val- 
ues of 1.04^Q;^g kpc and 6llg^ km s~^ for the distance and 
the velocity, as listed in Table 2. Figure 1 shows the resid- 
ual position offsets in each coordinate after subtracting 
the best fit proper motion, with sinusoids corresponding 
to the best fit parallax overplotted. It is apparent from 
the figure that the residual scatter in declination is signif- 
icantly worse than that in right ascension. The angular 
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separation between the calibrator and the pulsar is pri- 
marily in declination (Aacai = 0.1°, AS^ai = 3.3°), which 
may explain the excess scatter. Our estimates of position 
uncertainty do not currently account for the asymmetric 
calibrator throw, an issue we return to in §4. In tests of 
the fit after dropping one epoch at a time, the obtained 
parallaxes range from it = 0.87 mas to 1.02 mas, attesting 
to the robustness of the reported value of0.91±0.16 mas 
for the parallax. 

For B1929+10, too, a good fit was obtained with a re- 
duced — 1-06 for 7 degrees of freedom: the best-fit val- 
ues are listed in Table 1. The longer time baseline of the 
5 GHz data (2 years compared to 0.77 year for the 1.4 GHz 
data) results in a more accurate proper motion determina- 
tion, and yields a parallax tt = 2.76 ± 0.14 mas, which dif- 
fers from the Brisken et al. (2002) value of 3.02 ±0.09 mas 
at the 2(7 level. The covariances between {fXc/is) and tt 
are (—0.37, —0.01) for the 5 GHz fit, much lower values 
compared to (0.79,0.32) for the 1.4 GHz data. A joint 
fit to both data sets yields a parallax that is essentially 
unchanged from the 5 GHz value, but with smaller er- 
rors: TT = 2.77 ± 0.07 mas, with a reduced = 0.97 for 
15 degrees of freedom and covariances (0.15, 0.14) between 
(/Xq, jjsg) and TT. The zero-point offsets for the two data sets 
are expected to differ since they are referenced to different 
nodding calibrators at different frequencies (J1945-I-0952 
at 1.4 GHz and ICRF J192840.8-I-084848 at 5 GHz), but 
we obtain exceptionally good agreement. For the best joint 
fit, the difference in the zero-point offsets (Aa = 0.6 mas, 
Ai5 = 2.4 mas) implies that the absolute positions agree 
to ^ 2.5 mas. The residual position offsets after sub- 
tracting the best fit proper motion are shown in Figure 2, 
with the best fit parallax overplotted. The most probable 
value of the distance is 0.36llo;oo8 PC, and the velocity 
is 177^5 km s~^. Table 2 summarizes these and other 
derived parameters, and we defer a discussion of the ram- 
ifications of these measured parameters to §5. The joint 
fits reported here are robust compared to previous mea- 
surements reported for the same pulsar (Table 3), and the 
gradual convergence of the measurements is also gratify- 
ing. 

4. ESTIMATING ASTROMETRIC ACCURACY 

In phase-referenced observations, the phases of the tar- 
get visibilities are inferred from those of the nodding cal- 
ibrator. For each antenna, fluctuations in the ionosphere 
and troposphere are sampled on an angular scale equal to 
the separation between the target and calibrator sources. 
The level of fluctuation is characterized by the phase struc- 
ture function of the atmosphere (see, for example, Thomp- 
son, Moran & Swenson 2001, 13.1 and 13.3). For VLB ar- 
rays, the atmosphere is quite different at each antenna, and 
the uncalibrated differential phase for each visibility, which 
is a baseline-dependent quantity, cannot be interpreted as 
a simple excess in ionospheric path length or water va- 
por column density. It is reasonable, however, to expect 
that the effectiveness of phase referencing in calibrating 
target visibilities (and thus the quality of target images) 
depends on the angular separation between the target and 
the calibrator over which phase solutions have to be ex- 
trapolated. The temporal separation between scans on the 
calibrator and target introduces the time- variability of the 
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Fig. 1. — The parallax signature of PSR B0355-I-54 in right as- 
cension and declination, after subtracting the best-fit proper motion 
from the astrometric positions. Sinusoids corresponding to the best 
fit parallax tt = 0.91 mas are overplotted. 
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Fig. 2. The parallax signature of PSR B1929-I-10 in right as- 
cension and declination, after subtracting the best-fit proper motion 
from the astrometric positions measured at 5 GHz (closed circles; 
this work) and 1.4 GHz (open squares; Brisken et al. 2002). Si- 
nusoids corresponding to the best fit parallax tt = 2.77 mas are 
overplotted. 

atmosphere as another source of error, but this is modest 
in comparison, especially for short cycle times when phase 
connection is achieved across multiple scans. 

For the 5 GHz observations described here, B0329-I-54 
and B0355-I-54 shared the same nodding calibrator, and 
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Table 1 

Measured Positions, Proper Motions and Parallaxes 



Pulsar 


Q!2000^ 


^2000^ 


Ma 




TT 








(mas yr~^) 


(mas yr~-^) 


(mas) 


B0355+54 


03'^58'"53!71650 


+54°13'13'.'7273 


9.20 ±0.18 


8.17 ±0.39 


0.91 ±0.16 


B1929+10 


19'^32'"13!94970 


+ 10°59'32'.'4198 


94.03 ±0.14 


43.37 ±0.29 


2.76 ±0.14 


B1929+10 joint 


19^32'"13!94969 


+10°59'32'/4203 


94.09 ±0.11 


42.99 ±0.16 


2.77 ±0.07 



^Absolute positions extrapolated to epoch 2000.0 are accurate to 5 mas, which includes the 
uncertainties in the fit position and extrapolation based on proper motion (~ 1 mas), as well as 
absolute position uncertainties (< 3 mas) for the ICRF sources used in phase referencing. 



Table 2 
Derived Parameters 



Pulsar 


DM 

(pc cm~^) 


Distance 
(kpc) 


DnE2001 

(kpc) 


(km s~^) 


(rie) 
(cm~^) 


B0355+54 
B1929+10 


57.14 
3.176 




1.45 
0.33 


mtt 


0.052 
0.0087 



Table 3 

Measured Parallaxes for B1929+10 



Band 


TT 


Reference 


(GHz) 


(mas) 




0.4 


21.5 ±0.3 


1 


2.7 


< 4 


2 


2.2 


5.0 ± 1.5 


3 


1.4 


3.02 ± 0.09 


4 


5 


2.76 ±0.14 


This work 


1.4, 5 


2.77 ±0.07 


Joint fit 



References. — (1) Salter, Lync 
& Anderson (1979); (2) Backer 
& Sramck (1982); (3) Campbell 
(1995); (4) Brisken et al. (2002) 
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scans on the two pulsars were interleaved in time. But 
as described in §3, the quality of astronictry differed 
markedly for the two, leading to no useful result for one 
and a good sub-mas parallax for the other. The primary 
source of error is the 5.3° calibrator throw for B0329+54, 
as compared to only 3.3° for B0355+54. The observed 
position for the calibrator source J1934+1043 also var- 
ied when phase referenced to the nodding calibrator 2.4° 
away, while the positions for B1929-I-10 were well mod- 
eled by a proper motion and parallax after referencing to 
J1934+1043, 0.6° away. We quantify the random astro- 
metric error for each pulsar at each epoch by subtracting 
the best fit proper motion and parallax model and combin- 
ing the residuals for each coordinate in quadrature. The 
median error for each of the three target pulsars is plotted 
against calibrator throw in Figure 3, along with the an- 
gular position scatter for J1934-I-1043 about the median 
position. In each case, the central 50% range of the esti- 
mates is indicated with error bars. We note that this plot 
includes just 7 epochs, observed at only two declinations 
and three sun angles, with potentially unusual ionospheric 
activity at some epochs due to solar maximum. However, 
the trend of increasing median errors with calibrator throw 
is immediately apparent for the 5 GHz data. 
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Fig. 3. — Astrometric position errors as a function of the angu- 
lar separation between the target and the calibrator source. The 
points indicate the median scatter in the position of one calibra- 
tor measured with respect to another (outlined), or the offset of 
an astrometric target position from that predicted by the best-fit 
proper motion and parallax model. The error bars indicate the in- 
ner 50% of the range in the scatter, between the first and third 
quartilc. The 5 GHz data (open circles) are from this work. The 
1.6 GHz data (solid squares) include results from in-beam calibra- 
tion for PSR B0919-I-06 (left-most point; Chatterjce et al. 2001) 
and VLBI astroinetry of circumstcUar OH masers (Vlemmings et 
al. 2003). There arc between 4 and 12 samples at each calibrator 
throw, limiting the statistical significance of this plot at present, 
but the increase in astrometric error with calibrator throw is an 
unmistakable trend. Astrometric accuracy in phase referenced ob- 
servations is also better for 5 GHz data compared to 1.6 GHz data. 

For comparison, we also show astrometric data at 



1.4-1.6 GHz from prior projects. The parallax of 
PSR B0919+06 has been measured using in-bcam calibra- 
tion (Chatterjee et al. 2001), effectively using a calibrator 
throw of only 0.2°. The median residual random error, 
measured as above, is 0.5 mas. The parallaxes of cir- 
cumstcUar masers have been measured by Vlemmings et 
al. (2003) using standard phase referencing at 1.6 GHz, 
and the median residual errors range between 1-3 mas. 
The median errors are indicated in Figure 3 with solid 
squares, along with error bars showing the inner 50% 
range. Additionally, Vlemmings et al. (2003) performed 
relative astrometry between calibrator pairs, and we indi- 
cate the measured median relative scatter with outlined 
squares. Caveats similar to those listed for the 5 GHz 
data also apply here, since there are only between 4 and 
12 epochs used for each measurement. As in the case of 
the 5 GHz data, though, the astrometric error at 1.6 GHz 
generally increases with increasing calibrator throw. The 
errors are also significantly higher than those measured at 
5 GHz for comparable calibrator throws. 

Figure 3 empirically demonstrates the potential superi- 
ority of ordinary phase-referenced astrometry at 5 GHz, 
and emphasizes the need for additional techniques such as 
in-beam calibration (Fomalont et al. 1999; Chatterjee et 
al. 2001) or wide-band ionospheric calibration (Brisken et 
al. 2000) to obtain higher astrometric accuracy at lower 
frequencies. In §3, we noted that the position errors 
for B0355+54 in each dimension (Act, AS) appeared to 
be correlated to the calibrator throw in that dimension 
(Attcai, A(5cai)- This trend is observed for all the objects 
in Figure 3, at both 1.4 and 5 GHz. While the errors 
(Aa, A6) are unlikely to be completely independent of 
each other, the calibrator throw in each dimension is ap- 
parently an important predictor of the accuracy attain- 
able in that dimension. Large programs now underway 
will provide data on astrometric precision at currently un- 
sampled calibrator throws, as well as sampling a range of 
sun angles and a larger span of the solar cycle, enhanc- 
ing the utility of future versions of this plot both as an 
observation planning tool and as a probe of the tropo- 
spheric and ionospheric structure functions. Meanwhile, 
the current results reiterate the importance of using the 
smallest possible calibrator throw in phase-referenced ob- 
servations. Additionally, based on these and other recent 
results (Chatterjee et al. 2001; Brisken et al. 2002), we 
predict that parallaxes and proper motions should be mea- 
surable for the majority of pulsars within 2 kpc, and even 
for a some out to ~ 3 kpc, using either observations at 
5 GHz, or lower frequency observations with additional 
calibration techniques. Beyond this distance range, the 
required precision may run up against the limits imposed 
by potential calibrator source structure and variability at 
the 100 jias level. 

5. IMPLICATIONS: B0355-I-54 AND B1929-I-10 

The measured parallaxes and proper motions for 

B1929+10 and B0355+54 lead to estimates of the dis- 
tance (D) and velocity perpendicular to the line of sight 
{V±) for each pulsar. The most probable vahies for D and 
V± and the most compact 68% probability intervals are 
summarized in Table 2. The^istance estimated from the 
dispersion measure (DM = ne{s) ds) using the NE2001 
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model (Cordes & Lazio 2002) is listed for comparison. This 
model already includes a distance constraint for B1929+10 
based on the Brisken et al. parallax and overestimates the 
measured distance to B0355+54 by ~ 40%. Table 2 also 
lists the average electron density integrated over the line 
of sight, which is given simply by DM/D. The large dif- 
ference between the two lines of sight reflects the complex 
structure in the local interstellar medium, as discussed, 
for example, by Toscano et al. (1999) and Cordes & Lazio 
(2002). 

The structure of the local interstellar medium also af- 
fects the scattering properties of the pulsar signal, and 
measured scattering parameters can be used to constrain 
the distribution of scattering material along the line of 
sight (e.g. Chatterjee et al. 2001). Preliminary investiga- 
tions suggest that the scattering properties of B1929-I-10 
are affected by interstellar material which may be located 
quite close to the solar system. However, the velocity of 
the Earth is not negligible compared to the measured ve- 
locities of B1929+10 and B0355+54; and we postpone a 
detailed analysis of the interstellar scattering properties of 
these pulsars to future work. 

The astrometric results discussed here also lead us to 
re-examine, in turn, the birth site for B1929+10 based on 
its kinematics, the nature; of its observed X-ray emission, 
and constraints on the cooling curve for each of the two 
neutron stars. 

5.1. The Birth Site of B1929+10 

Based on numerical simulations, Hoogerwerf, de Bruijne 

& de Zeeuw (2001) have suggested that PSR B1929+10 
(which they refer to as J1932+1059) and the nmaway O- 
star ( Ophiuchi were part of a binary system, which be- 
came unbound when the pulsar progenitor exploded in a 
supernova ~ 1 Myr ago in the Upper Scorpius subgroup 
of the Sco 0B2 association. They attribute the kick ve- 
locity of the pulsar and the large space velocity of ( Oph 
20 km s""'^) to the same supernova- induced binary dis- 
ruption event. As seen in Figures 4 and 5 of Hoogerwerf, 
de Bruijne & de Zeeuw (2001), only a limited region in 
the parameter space of pulsar parallax, proper motion and 
radial velocity allows the pulsar and ^ Oph to approach 
within 10 pc of each other and 10 pc of the center of the 
Upper Scorpius association. For these successful simula- 
tions (0.14% of the total), the implied kinematic age of 
B1929-I-10 is ~ 1 Myr, comparable to its spindown age of 
3 Myr. 

The distance and velocity measurements underpinning 
such a birth scenario are based on a "dearth of data" 
(Campbell 1995), and the new values for the proper motion 
and parallax of B1929+10 that we report here make the 
scenario extremely unlikely. For tt = 2.76 mas, a high ra- 
dial velocity (> 300 km s~^) is required, while our proper 
motion measurements arc inconsistent with the parame- 
ters used in the successful simulations by > lOcr. There- 
fore, B1929-I-10 is probably not the former binary com- 
panion of ( Oph. We note, however, that the origin of the 
pulsar in the Upper Scorpius region is not excluded, and 
the pulsar could still have had its birth in a supernova in 
this region 1-2 Myr ago (de Geus 1992). The isolated neu- 
tron star RX J1856.5— 3754 had initially been proposed as 
a potential binary companion to ( Oph in the Upper Scor- 
pius association (Walter 2001), but the revised distance to 



the neutron star (Kaplan, van Kerkwijk & Anderson 2002; 
Walter & Lattimer 2002) rules this scenario out, since the 
closest approach by RX J1856.5— 3754 to the association 
center occurred only ~ 0.5 Myr ago. Thus, for ( Oph, the 
alternative birth scenario in the Upper Centaurus Lupus 
subgroup 2-3 Myr ago (van Rensbergen, Vanbeveren, & 
de Loore 1996) looks more plausible. 

5.2. The Nature of X-ray Emission from B 1929+ 10 

X-ray emission from PSR B1929+10 was detected in ob- 
servations with the Einstein (Helfand 1983), ROSAT{Yan- 
copoulos, Hamilton & Helfand 1994), and ASCA (Wang & 
Halpern 1997) observatories. The latter two datasets allow 
for spectral analysis that can be used to constrain proper- 
ties of the X-ray emitting region. Yancopoulos, Hamilton 
& Helfand (1994) fit the ROSATdata to a model for black- 
body emission, deriving a temperature of roughly 3 x 10^ K 
(hereafter, 1 MK = 10^ K). From ^S'C^-SIS data, Wang & 
Halpern (1997) find comparably successful fits to a single 
blackbody (T = 5.14 ±0.53 MK) and to a power-law spec- 
trum. Finally, Wozna, Kuiper & Hermsen (2003) perform 
spectral fits to combined RO SAT and ASCA data (using 
in this case the ASCA-GIS detector) spanning the energy 
range 0.1-10 keV; they find acceptable fits to model spec- 
tra consisting of either a single power law or two black- 
body components with temperatures of 2 and ^ 7 MK. 
New observations will be needed to ascertain whether the 
emission is thermal or nonthermal in nature. If thermal, 
the existing data alone would rule out an origin for the 
X-rays in "cooling" emission from the entire surface of 
the neutron star, both because the implied temperature 
is unexpectedly high for such an old pulsar (characteristic 
age ~ 3 Myr) and because the emitting area would be un- 
physically small. Indeed, as discussed by Wang & Halpern 
(1997), the inferred emitting area is nearly two orders of 
magnitude smaller than theoretical predictions even for 
the size of a heated polar cap assuming a purely dipolar, 
star-centered magnetic field. Wang, Ruderman, Halpern 
& Zhu (1998) elaborate on this finding, arguing that a 
small observed cap size may be evidence for a low heating 
current from the magnetosphere, or an off-center dipolar 
geometry, perhaps as a result of crustal plate motion on 
the neutron star surface. 

The implications of a new distance determination de- 
pend on whether the emission from the neutron star is 
thermal or nonthermal in nature. If the emission is as- 
sumed to be thermal, and accurately represented by the 
single-blackbody result of Wang & Halpern (1997), our 
new distance measurement implies a bolometric luminos- 
ityi L = (2.74 ± 0.37) x 10^° ergs s"\ and an emitting 
area A = (6.8 ± 3.0) x 10^ cm^, where formal la statisti- 
cal uncertainties in spectral fitting (i.e., neglecting possible 
systematic effects due to calibration uncertainties) and dis- 
tance have been propagated, and isotropic emission over 
47r sr has been assumed. The new parallactic distance 
does not, by itself, alleviate the apparent discrepancy be- 
tween observed and expected polar cap sizes, increasing 
the inferred emitting area by a factor (361/250)^ ~ 2. 

Alternatively, the thermal fits of Wozna, Kuiper & 

^ Unless otherwise noted, all luminosities, temperatures, and emit- 
ting areas are quoted for a distant observer. Where adjustments 
are made for rcdshift at the neutron star surface, canonical values 
M = 1.4 Mq and R = 10 km, are assumed. 
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Hcrmscn (2003), to ROSAT and ASCA data simulta- 
neously, imply omission luminosities of Li = (1.85 ± 
0.30) X 103° ergs s'^ and L2 = (1.72 ± 0.30) x lO^^ 
ergs s~S with areas A-^ = (1.88 ± 0.35) x 10^ cm^ and 
A2 = (1.34 ± 0.36) X 10^ cm2, for the 2 MK and 7 MK 
blackbody components, respectively. If this spectral model 
holds, the intrinsic emitting area of the low- temperature 
component, 

Ai"* = Ai(l - 2GM/Rc^) = 1.1 X 10^ cm^, (1) 

comes within a factor of three of the canonical polar-cap 
size defined by the last open field lines in the pulsar mag- 
netosphere, 

Ape = = 2.9 X 10^ cm\ (2) 

where P is the pulse period, an intriguing similarity given 
uncertainties in magnetospheric geometry, the beaming 
solid angle, and likely distortion of the surface blackbody 
spectrum due to propagation of the radiation through a 
stellar atmosphere. 

If, instead, the observed X-ray emission is interpreted 
correctly by the nonthermal model of Wozna, Kuiper & 
Hermsen (2003), then the 0.1-10 keV luminosity implied 
by our new distance measurement is _LpL = (9.6 ± 0.8) x 
10^° ergs s~^; in the 2-10 keV band, the same model im- 
plies a luminosity 3.1 x lO"'^ of the pulsar's spin-down 
power, comparable to the nonthermal X-ray emissions 
of other pulsars (Possenti et al. 2002). 

Finally, we address the implications of our new distance 
determination on stellar cooling models for B1929+10. 
Following the procedure adopted by Yancopoulos, Hamil- 
ton & Helfand (1994), we derive upper bounds on the sur- 
face temperature by requiring that whole-surface "cool- 
ing" emission not exceed the observed soft-band (0.09- 
0.30 keV) ROSAT countrate. The result is sensitive to 
the assumed absorption column, njj. For the parallac- 
tic distance and the smallest likely nn, 1-5 x 10^° cm~^, 
a temperature T < 0.24 MK is required, corresponding 
to a luminosity limit L < 4 x 10"^° ergs s~^, roughly an 
order of magnitude larger than typical theoretical expec- 
tations (~ 5 X 10^^ ergs s~^) for residual heat retained 
since the formation of such an old pulsar (Tsuruta et al. 
2002). If uh is allowed to be as large as 6.5 x 10^° cm~^, 
the 1(7 upper limit for thermal models derived by Wozna, 
Kuiper & Hermsen (2003), an even hotter surface is al- 
lowed: T < 0.35 MK and L < 2 X 10^^ ergs 8"^ We 
therefore conclude that available soft X-ray measurements 
do not usefully constrain leading cooling models. 

5.3. High Energy Emission from B0355+54 

PSR B03554-54 was detected by the ROSAT observa- 
tory (Slane 1994). Too few photons were acquired to 
permit detailed spectral modeling, but the gross distri- 
bution of photon energies suggested nonthermal emission. 
A recent observation with the XMM-Newton telescope has 
confirmed this conclusion (Kennea et al., in preparation). 
Adopting a distance of 2.1 kpc, Slane (1994) inferred a 
nonthermal luminosity (0.1- 2.4 keV) of 1.0 x 10^^ ergs s"-'^; 
scaling to the parallactic distance of 1.04 kpc results in 
a revised luminosity of 2.4 x 10^^ ergs s~^. The avail- 
able XMM-Newton data also provide a Icr upper bound 
on the surface cooling luminosity of the neutron star, 



L < 4.5 X 10'^^ ergs s^^ (J. Kennea, private communi- 
cation). As with B1929+10, this limit is consistent with 
leading cooling models (Tsuruta et al. 2002). 

The new distance measurement for B0355-I-54 has simi- 
lar implications for studies of the 7-ray emission from the 
pulsar. The EGRET instrument aboard CGRO set an 
upper limit to the flux of 0.1-10 GeV photons that im- 
plies a conversion efficiency r]j = Lj/E < 0.2 (Nel et al. 
1996), if the 7-ray beam crosses the line of sight. Our par- 
allactic distance, nearly a factor of two smaller than the 
dispersion-derived distance, further reduces the apparent 
7-ray efficiency, rj-y < 0.05, placing B0355-I-54 at odds with 
the observed trend of increasing 77^ with spindown age. 
For example, PSR B1055— 52, with a characteristic age 
similar to B0355+54, has rj^ ~ 0.7. If the trend of 7-ray 
efficiency increasing with age holds, the likely conclusion 
is that the 7-ray beam for B0355-I-54 does not cross our 
line of sight, as predicted in the outer-gap beaming model 
of Yadigaroglu & Romani (1995). 

6. CONCLUSIONS 

We have presented results from VLBA astrometric ob- 
servations of three pulsars at 5 GHz. For B0329-I-54, the 
systematic errors are large enough that we obtain only an 
upper limit to the parallax, tt < 1.5 mas, and proper mo- 
tions consistent with the Brisken et al. (2002) value. The 
large angular separation between the calibrator and pulsar 
(5.3°) is the most probable cause for the large systematic 
errors. With the same calibrator and at the same observa- 
tion epochs, but with a smaller angular throw (3.3°), we 
obtain a proper motion and parallax for B0355+54, with 
an implied D = 1.04^°;^^ kpc and V± = 61^^^ km s"!. 

For B1929-I-10, we demonstrate the consistency of our 
absolute astrometry with that of Brisken et al. (2002) to 
within 2.5 mas. The longer time baseline of the 5 GHz 
data permits a better proper motion measurement while 
reducing the covariance between fi and tt. A joint fit to 
both data sets produces a parallax identical to that mea- 
sured with the 5 GHz data alone (but with smaller un- 
certainties), implying a distance D = 0.36llo;oo8 ^P*^ 
velocity V± = 177ti km s"^. The results for B1929-M0 
illustrate the need for both a time baseline long enough 
to minimize the covariance between n, tt, and a sufficient 
number of epochs so that the parallax signature is well 
sampled. 

Based on the observed astrometric errors in our obser- 
vations, as well as the astrometric uncertainties measured 
by Chatterjee et al. (2001) and Vlemmings et al. (2003) 
at 1.4 GHz, we quantify the astrometric error as a func- 
tion of calibrator throw and observing frequency. The as- 
trometric precision is correlated with the angular throw 
between the calibrator and the target, with lower errors 
at 5 GHz compared to 1.4 GHz. Figure 3 provides a first 
attempt at quantifying this dependence on an empirical 
basis, although we note that the random nature of these 
errors implies that a value of parallax smaller than the 
predicted astrometric error can be still be measured given 
enough measurement epochs. The calibrator throw in each 
dimension may serve as an important predictor of the ac- 
curacy attainable in that dimension. Future observations 
will allow us to confirm or refute this apparent correlation. 

Finally, we have investigated the implications of the 
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new distance and transverse velocity determinations for 
B1929+10 and B0355+54. We conclude that B1929+10 
and the runaway 0-star Ophiuchi do not approach 
within 10 pc of each other when their paths are traced 
back in time, given the revised distance and velocity for 
the pulsar. Thus, contrary to the suggestion by Hooger- 
werf, de Bruijne & de Zeeuw (2001), they arc unlikely to 
have been members of a binary system that was disrupted 
by the birth supernova of the pulsar. A review of X-ray 
data allows us to constrain the luminosity of the two pul- 
sars in soft X-rays. The limits on the thermal emission 
from both objects are consistent with (and unconstraining 
for) leading neutron star cooling models (Tsuruta et al. 
2002). If the X-ray emission from B1929-hl0 is thermal 
in nature, the emitting area for X-rays is roughly consis- 
tent with the canonical size for a heated polar cap. For 
B0355-I-54, the inferred conversion of spin-down power to 
gamma-ray luminosity is low, which may imply that its 
7-ray beam does not intersect our line of sight. 

The implications that wc discuss here illustrate the 
range and diversity of the applications of new distance 
and velocity measurements for neutron stars. Most pul- 
sars are weak at 5 GHz, but we expect that a few others 
will be observable with the VLBA alone. The addition of 
telescopes with large collecting areas such as Arccibo and 
the Green Bank Telescope to the array will allow us to 
further enlarge the sample, and explore the consequences 
at multiple wavelengths. 
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